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Abstract 

The initial bifurcations in rotating Rayleigh-Benard convection are 
studied in the range of dimensionless rotation rate < Q < 2150 for an 
aspect-ratio-2.5 cylindrical cell. We used simultaneous optical shadow- 
graph, heat transport and local temperature measurements to deter- 
mine the stability and characteristics of the azimuthally-periodic wall 
convection state. We also show that the second transition corresponds 
to the onset of bulk convection. Our results for critical Rayleigh num- 
bers, precession frequencies and critical mode numbers agree well with 
theoretical results. The dynamics of the wall convection state can be 
described by a complex Ginzburg-Landau amplitude equation. 

1 Introduction 

Studying the effects of external rotation on thermal convection has attracted 
significant experimental and theoretical interest [p], 0, |3], |J. Because of its 
general occurrence in geophysical and oceanic flows, it is important to under- 
stand how the Coriolis force influences the structure and transport properties 
of thermal convection. Rotating thermal convection also provides a system 
in which to study hydrodynamic instabilities 0, [5], || , pattern formation and 
spatio-temporal chaos in nonlinear dynamical systems @, |7j, and transitions 
to turbulence || [|. 



Rayleigh-Benard convection consists of a fluid layer confined between two 
horizontal plates and heated from below. External rotation is imposed by 
rotating the convection cell about a vertical axis. Most theoretical studies 
of this system have assumed a laterally infinite cell geometry Jl|. In experi- 
ments, however, there are always side walls. Moreover, to achieve sufficient 
laboratory control over the experimental parameters and obtain large dimen- 
sionless rotation rates, the horizontal dimensions of the convection cells are 
usually comparable to their vertical depths. This raises the question whether 
the theoretical predictions for infinitely extended systems are relevant for real 
experiments. It has indeed been found that the presence of lateral bound- 
aries change significantly some aspects of thermal convection in the presence 
of rotation^], f|, 0. In this paper, we present experimental results that il- 
lustrate clearly the influence of the lateral boundaries. We confine ourselves 
to cylindrical convection cells and characterize the geometry by the aspect 
ratio which is the ratio of the radius to the depth. 

There have been several experimental and theoretical investigations on 
rotating Rayleigh-Benard convection. One of the early quantitative exper- 
iments by Rossby|| showed significant deviation from the predictions of 
theories for a laterally infinite system. His measurement of heat transport 
showed a slow increase over the pure conduction state at temperature differ- 
ences much smaller than the expected onset of convection for an infinitely 
extended system. There have been many explanations offered to remedy 
this discrepancy. The first one that correctly attributed the increase in ther- 
mal transport to a side wall convection state was proposed by Buell and 
Catton||. They investigated the stability of the conduction state in a lat- 
erally finite cylinder in the presence of rotation. Contrary to conventional 
wisdom that side walls suppress thermal convection, they found that for a 
large enough rotation rate, some side wall convection states have a critical 
control parameter that is much lower than the value for the laterally infi- 
nite case. These non-axisymmetric wall states are azimuthally periodic and 
have maximum amplitude near the side wall. Zhong et a/Q later reported 
experimental confirmation of the wall convection states in a unit aspect ratio 
cell. They observed, however, that these asymmetric states precessed slowly 
counter to the direction of rotation. In addition to the first bifurcation to 
the wall states, there was a second bifurcation when the central region be- 
came filled with vortex-like structures, accompanied by an increase in the 
thermal transport and in the thermal noise measured in the bottom plate. 
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Goldstein et a/[|Hj] carried out a linear stability analysis which allowed for 
the experimentally observed precession and found good agreement with the 
experimental results. Ecke et al\Tl\ explained the traveling wall states by 
a broken azimuthal reflection symmetry and predicted that the precession 
frequency would go linearly to zero as the rotation rate went to zero|TT] . 



Our experiment, using an aspect ratio 2.5 cell, was designed to explore 
the aspect ratio dependence of the initial bifurcation from the conduction 
state to the wall convection state. We chose to increase the aspect ratio only 
moderately relative to the previous cell of Zhong et al @] because the range 
of accessible dimensionless rotation rates (limited by the condition that cen- 
trifugal effects be small) for the larger aspect ratio cell overlapped a large 
fraction of the parameter range used in the previous work. Furthermore, 
although the physical aspect ratio is considered small for nonrotating con- 
vection, rotation increases the effective aspect ratio significantly. In addition, 
since the side-wall state forms around the azimuth, the number of structures 
that fit in that geometry is roughly the ratio of the perimeter to twice the 
depth. For F = 2.5 this ratio is 7.9. Thus, provided the dimensionless ro- 
tation rate is sufficiently large, one can consider the conditions to be in the 
moderate to large aspect ratio limit and our measurements presented here 
together with the previous work for T = 1 suffice to characterize the aspect 
ratio dependence of the observed bifurcations. We discuss this in more detail 
later. 

In this work we present measurements of the critical Rayleigh numbers 
and critical wavenumbers for the initial transition to a side-wall traveling 
wave state and the second transition to bulk convection. We find that the 
traveling wave state can be described quantitatively by a one- dimensional 
complex Ginzburg-Landau (CGL) equation. In addition we find that the 
traveling wave state does not depend on the overall geometry of the convec- 
tion cell. Rather it exists on any wall regardless of the magnitude or sign of 
the wall curvature. Further, the conductivity of the wall is shown to strongly 
influence the critical Rayleigh number for the bifurcation to traveling waves. 
Finally, the second transition to bulk convection is shown to occur as a spa- 
tially seperate transition owing to the negligible amplitude of the wall state 
in the center of the cell. 
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2 Experiment 



The apparatus used in these convection studies was descibed in detail by 
Zhong et a/Q. The bottom plate of the convection cell was a nickel-plated, 
mirror-polished, copper plate with a thermofoil heater attached to the back. 
The top plate was a 1/8 inch thick, optical quality sapphire window. The 
side wall was made of 0.31 cm thick plexiglasss, defining an active circular 
area with a radius r = 5.00 cm. The depth d of the cell was 2.00 cm, giv- 
ing an aspect ratio T = r/d = 2.5. We used an afocal optical shadowgraph 
technique to visualize the convective flow field. Because of the depth of the 
cell (hence a very small temperature gradient), the shadowgraph method 
was not sensitive enough to observe the flow near the first bifurcation, espe- 
cially when the rotation rate was small. For local temperature measurement, 
two thermistors were embedded at half-height positions in the side wall with 
known angular separation. Since the wall states traveled along the circum- 
ference, the thermistors measured the temperature field sweeping by. From 
these data, the amplitude, frequency and mode number of the traveling wave 
were extracted. The heating was provided by applying a constant voltage 
to the bottom heater, and a circulating, temperature-regulated water bath 
outside the top sapphire window provided cooling and constant top temper- 
ature to better than 0.5 mK rms. Since the convection cell rotates relative 
to the regulated water bath, the top-plate temperature is azimuthally very 
uniform on the time scale relevant to the dynamics of the system. Ther- 
mistors embedded in the bottom plate and near the top plate measured the 
respective temperatures to a precision better than 0.1 mK. Rotation was 
provided with a stepper motor with rotation frequency constant to ±0.1%. 
Because the shadowgraph optics was at rest in the laboratory, a shaft en- 
coder strobed the video camera at each revolution of the cell. The image 
digitized via a digital frame grabber was thus synchronized with the rotating 
cell. Since the dynamics of the convection patterns were slow compared to 
the rotation period, we did not lose any detailed information by recording 
the images once per revolution. 

The control parameters of the experiment are Rayleigh number R and 
dimensionless rotation rate Q, defined as R = agd 3 AT/Ku and Q = Q^d 2 /u, 
where a, k and v are the thermal expansion coefficient, thermal diffusivity 
and kinematic viscosity of the working fluid, g is the gravitational acceler- 
ation, AT is the bottom-top temperature difference and Qd is the angular 
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rotation frequency. Another parameter that characterizes the fluid is the 
Prandtl number P = v/k, which is about 6.4 for our choice of water at a 
mean temperature of about 23.5°C. A reduced bifurcation parameter is de- 
fined as e = (R — R c (fl)) / R c (fl) , where R c (£l) is the critical Rayleigh number 
for the onset of the first convection state at a particular Q. As we will dicuss 
below, the first bifurcation can be to different states depending on Q. There- 
fore, we denote the critical Rayleigh number for the side-wall state as R cs 
and that for the onset of bulk convection as R c b. The effect of the centrifu- 
gal force can be characterized by the ratio of the centrifugal-to-gravitational 
force, Q 2 D r/g, which is less than 0.12 for our maximum rotation rate Qd ~ 
4.9 rad/sec (ft « 2145). 

The experimental protocol we used was as follows: Qd was held con- 
stant for the desired f2 at AT C ; the heating was quasistatically ramped 
to each measurement point. After several vertical thermal diffusion times 
t k — d 2 1 k ~ 2760 sec, measurement of temperatures were made spanning 
at least one r K . Shadowgraph images were recorded at fixed time intervals 
during the temperature measurements. We followed the dynamics at sev- 
eral parameter values for longer than 40r K , and no change of the convection 
states was found for the parameter range studied. There is one exception 
concerning the stability of the wall state which will be discussed later. 

For determining the onset of the wall convection state, two different mea- 
surements were used for large enough values of Q so that the transition to the 
wall state was the first bifurcation. For smaller f2 only the second method 
was used. The first was the measurement of heat transport, expressed in Nus- 
selt number Nu which is the ratio of the heat transported by the convection 
state to the heat transported by conduction alone. In the conduction state 
Nu = 1 whereas just above the onset of convection, Nu increases linearly 
with increasing R. The intercept of the linear section of Nu versus R with 
Nu = 1 determines R cs . The second method used the temperature measure- 
ments from the side wall thermistors T m( i\ and T m d 2 - By fitting the oscillatory 
signals to a sine function, we extracted the amplitude 5T, which rises above 
the onset of the side wall state as (R — /^(ft)) 1 / 2 . The two methods yielded 
values of R c that agreed to within experimental uncertainty, usually better 
than 0.2%. The second method also gave us the frequency and mode number 
of the traveling wave. 

The second bifurcation was determined by the Nusselt number measure- 
ment. We did, however, find two other characteristics of the bulk convection 
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state that could be used to determine the onset, namely, the thermal noise 
in the bottom temperature and the amplitude of the shadowgraph images. 
These results will be discussed in detail later. 

3 Experimental Results and Comparison with 
Theories 

We present our experimental results on an aspect ratio 2.5 cell. When ap- 
propriate for showing variation as a function of aspect ratio, we also include 
results for the unit aspect ratio cell from Zhong et aZQ. Comparison with 
theories are made whenever theoretical results are available. 

3.1 Asymmetric wall convection states 

As the temperature difference between the bottom and top plates increases 
at fixed Q (Q actually changes slightly in the experiment since we held 
constant; typically, Q varies less than 0.4% for the range of Rayleigh num- 
bers we used), the first bifurcation from the quiescent conduction state is 
to an azimuthally periodic wall convection state for moderate to high rota- 
tion rates. This convection state precesses uniformly in the rotating frame, 
counter to the direction of external rotation and with maximum amplitude 
near the side wall. The amplitude grows above the onset as e 1//2 and the 
frequency lo k changes linearly with e and has a finite value at onset. This 
indicates that the bifurcation is a forward Hopf bifurcation. The unique di- 
rection of precession arises from the weak breaking of azimuthal reflection 
symmetry by the external rotation (Ecke et q/)||ll||. The e dependence of 
the amplitude and the frequency of the wall convection state is adequately 
characterized in Zhong et al[Q . We concentrate on the Q dependence of the 
critical frequencies and mode numbers, and in particular, on their behaviors 
as Q goes to zero or becomes large. We also discuss the linear stabilities for 
different mode-number wall convection states. 

We have measured the critical Rayleigh number of the wall convection 
states in a wide range of Q. The results are displayed in Figure [l]. The 
dashed line is from the calculations of Goldstein et al for perfectly insulating 
side walls [|10|, corresponding reasonably well with the experimental boundary 
conditions (plexiglass side wall). The agreement between experiments and 
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Figure 1: Parameter-space diagram of rotating Rayleigh-Benard convection 
in small aspect ratio cylindrical containers. Measured bifurcation points 
from conduction to wall convection (A) and bulk convection states (A) for an 
aspect ratio 2.5 cell and corresponding bifurcation points (• and o) measured 
in a unit aspect ratio cell|| are shown. Rossby's results (+) for aspect ratios 
1.4-50, Chandrasekhar's linear stability calculation for a laterally infinite 

system ( — ), and Goldstein et al's result ( ) for a unit aspect ratio cell 

with insulating side walls |[lC|| are also shown for comparison. 

theory is very satisfactory. The results of Zhong et al for T = 1 are also 
shown in Figure |l|. For completeness, we have included experimental data of 
Rossby|J using cylindrical convection cells with a range of aspect ratios from 
1.4 to 5.0. There is not much variation in R cs among these cell geometries for 
Q > 100. We believe that the wall-state bifurcation curve is quite insensitive 
to further increases in T. The other bifurcation line and the cross-over of the 
two bifurcation lines in {R, Q} parameter space will be discussed later. 

Figure shows the measured precession frequency at onset u k0 scaled by 
t k for the wall states versus Q in a semi-log plot. The inset shows the linear 
dependence of the critical frequency on Q for Q < 150. The large Q behavior 
of uj kQ suggests a nearly logarithmic dependence on Q. Notice that for high 
fls, the frequencies for aspect ratios 1 and 2.5 cells approach each other. 
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Figure 2: Comparison of experimental and theoretical results on the preces- 
sion frequencies of the wall convection states. Data for aspect ratio 2.5 (o) 
and 1 (•) cells, Goldstein et aPs results ( — ) for a unit aspect ratio cell, and 

the results of Kuo and Cross ( ) for a planar wall state. The inset is a 

linear plot of the frequencies at small rotation rate and the dashed line is a 
fit to the data. 
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Table V. Values of Q, R CS} uj k0} m c , A c , R^, and k c f,. 



The solid line in the figure comes from a recent calculation by Goldstein et 
al for a unit aspect ratio geometry |]3L 0f| . Their calculation also shows that 
for moderately high Q, the change in frequency with increasing aspect ratio 
is small for T > 1. The dash-dotted curve is the result for infinite aspect 
ratio cells with insulating rigid side wall and free-free top-bottom boundary 



conditions from Kuo and Cross 12 



The critical mode numbers m c of the wall states were measured by qua- 
sistatically increasing AT from below. To scale out the aspect ratio depen- 
dence and identify the trend as both T and Q increase, we convert the mode 
numbers to dimensionless azimuthal wavelengths scaled by the depth d, so 
that A c = 27iT/m c . In Table [I] we show data for the onset of the side- wall 
state denoted R cs , u K o, m c , A c , and the critical Rayleigh number R c b and 
the critical wavenumber k c j, for bulk convection for the full range of Q values 
we explored for the T=2.5 convection cell. Figure |3| shows A c versus Q for 
T of 1 and 2.5. The numbers in parentheseses near the data points are the 
corresponding mode numbers. The stepped solid line is from Goldstein et 
a/[[nj for T = 1. The large Q behavior of A c suggests that it approaches a 
value near 1 or, in other words, that the critical mode number approaches a 
finite constant. This implies that the depth of the convection cell is setting 
the spatial scale of this wall state. 

We can understand the behavior of m c (Q) in the following manner: using 
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Figure 3: The critical wavelength for wall convection states with T=2.5 (o) 
and T=l (•). The bracketed numbers are the corresponding mode numbers. 
The solid line is from Goldstein et al for a unit aspect ratio cell and the 
dash-dotted curve is from Kuo and Cross for a planar wall state. 
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a physical argument by Manneville [ D^ ] , we know that to minimize the dissi- 
pation losses due to vertical and horizontal shear, the convection rolls should 
have similar width and height. In the case of wall convection in the presence 
of rotation, since 1) the Coriolis force only serves to skew the particle tra- 
jectories without providing kinetic energy, and 2) there is additional viscous 
damping in the vertical direction because the wall states lie close to the side 
wall, the convection structure should have a width comparable to but nar- 
rower than if it was not near a rigid wall. This leads to a critical wavelength 
in the streamline plane less than 2 for the wall states. In the bulk convection, 
external rotation gives rise to increasing apparent wavenumbers(measured as 
the distance between the periodic structures) as fl increases JL4|. In our ex- 
periment, since the particle trajectories near the side wall have to be nearly 
parallel to the wall, the measured wavelength is in the streamline plane. 
From the radial basis functions for the velocity and temperature perturba- 
tion field, which are mth order Bessel functions J m (kr), we can estimate the 
width of the convection structure in radial direction to be 5r « X m /k, where 
X m fa 3 is the width of an oscillation period for Bessel functions, and k has 
to satisfy a dispersion relation [|K|]. Since k increases with fl, the radial width 
of the wall states decreases as fl increases. We can reasonably argue that 
the wavelength at the middle of this width is 0(d), independent of rotation. 
Then m c {fl) = 2vr(r - 5r/2)/\ « 2vr(T - X m /kd), which gives for T=2.5 an 
asymptotic critical mode of m c « 16. This explains why the critical mode 
number increases with fl but saturates at about 2irT. 

In addition to the quantitative measurements presented above we per- 
formed a series of qualitative investigations to test the relevance of different 
cell geometries and side wall boundary conditions. In one arrangement a 
circular plexiglass wall of half the radius of the convection cell was place 
concentrically inside the cell, thus forming an annulus of width to depth ra- 
tio 1.2. Upon increasing R above R cs , wall states formed on all sides of the 
walls. In the annulus there was a traveling wave (TW) on the outer rim that 
traveled opposite to the direction of rotation while on the inner rim the TW 
traveled in the same direction as the rotation direction. In another arrange- 
ment, horizontally straight plexiglass walls were placed inside the cell, one 
side of which was covered with a thin copper tape providing an effectively 
conducting side wall boundary condition. TWs formed on the insulating side 
wall at R cs but were suppressed on the conducting wall. The direction of the 
TWs is fl x n where fl is the rotation vector and n the the normal vector 
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Figure 4: Marginal stability boundary in parameter space of R and azimuthal 
wavenumber k (•) and u^o vs k (q) for Q = 544 with T = 2.5. The solid line 
is a nonlinear third-order polynomial fit to the marginal stability data, the 
short-dashed line is the finite-size Eckhaus boundary, and the long-dashed 
line is a linear fit to the frequencies. 



of the wall into the fluid. A straight wall can thus be considered as the side 
wall for an effectively infinite convection cell. Kuo and Cross have recently 
calculated the properties of this straight wall state using free-free top and 
bottom boundary conditions and rigid boundary conditions at the lateral 
wall.|12| 

The amplitude of the wall convection state is concentrated near the side 
wall and behaves rather like a quasi-one- dimensional periodic pattern along 
the circumference. At Q = 544, we prepared several states with different 
mode numbers by jumping to a high AT and quickly dropping back to just 
above the critical AT. We then measured the amplitude of the traveling wall 
state and determined the marginal stability boundary by extrapolating the 
squared amplitude versus R to zero. Figure f| shows this marginal stability 
boundary and the corresponding precession frequency measured for Q = 544 
in the T = 2.5 cell. From fits to the data we obtain R c = 44217, k c = 4.65, and 
^kOc = 28.3. The actual switching of mode number toward the critical one 
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Figure 5: Temperature amplitude squared A 2 (•) and frequency u K (A) vs 
R for the transition from mode 13 to mode 12. 



occurred at finite amplitudes, suggesting an instability other than marginal 
instability. To identify this instability, we prepared an m = 13 state at R 
= 4.63 timeslO 4 and ramped down to R cs {m = 12) over about 80 hours 
making measurements at intervals of 5R ~ 0.002. We define e in terms of 
R c = 44217 for the marginal stability curve. At R = 4.547 x 10 4 , e = 0.0283 
with R c = 44217 for the marginal curve, traveling wave's amplitude began 
to deviate from a linear dependence on R (see Figure [5|). The transition to 
the critical mode number m c = 12 took place 27 hours (32r K ) later. The 
extrapolated marginal stability of the m = 13 mode is at R = 4.469 x 10 4 
(e = 0.0107). Fig. | also shows the measured precession frequencies scaled 



Because of the quasi one-dimensionality of the traveling wall states, it 
should be described by a 1-D complex Ginzburg-Landau (CGL) equation, 



Since the precession is opposite to the rotation direction frequencies are neg- 
ative. We have measured all the coefficients in the CGL equation at — 544, 



by t k . 




d 2 A 
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and they are c = 0.6, a = 0.6, -0.01 < c 3 < 0, r = 0.027, s = 2.1, £ = 
0.192, and go = 0.8, scaling time by d 2 / k and length by d. In this equation, 
e is defined relative to R c = 44217 and the modulation wavenumber rela- 
tive to k c =4.65, both obtained from the fit to the marginal stability curve, 
Fig. |j. Using these definitions we obtain the linear coefficients as follows: 
£o is the quadratic curvature of the marginal stability curve, Fig. [| ; r is 
obtained from transient measurements of the oscillation amplitude; So is the 
slope of the frequency at the critical wavenumber, Fig. |5|; the combination 
c\ — Co is determined from the curvature of the frequency dispersion, Fig. § 
the combination c 3 — cq comes from the variation of frequency with e. To 
resolve the remaining degree of freedom we measure the transient frequency 
response. To double check our result that c 3 is small we use a modulation 
technique in which we sinusoidally vary the bottom plate temperature with a 
period larger than 5.2r K . This method allows us to probe several parameters 
at once and to develop better statistics than is possible for the single tran- 
sient measurements. These results are consistent with the values obtained 
through the other static and dynamic measurements. The nonlinear coeffi- 
cient is obtained by normalizing the local probe temperature with the heat 
transport using Nu — 1 = \A\ 2 R/R C and with an area correction since the 
traveling- wave state is confined to a region near the wall of order the depth 
of the fluid (we take a radial width d so go = (dNu/de)(T 2 /(2T — 1)). We 
can compare our linear coefficients, co,ci,r ,so, and £0 with predictions of 
the linear theory for our particular experimental conditions Jl7| which yield 
c = 0.7, Ci = 0.8, r = 0.022, s = 1.8, and £ = 0.21, in reasonable agree- 
ment with the experimental values. As in the experiments the CGL equation 
is not exact so there are finite size uncertanties in the numerical coefficients. 



Infinite-system, nonlinear calculations |12| for a traveling-wave convection 
state on a flat wall for f2=540 and with free-free top-bottom boundary con- 
ditions yield c = 0.8, c x = 0.45, c 3 = 0.32, r = 0.022, s = 2.1, f = 0.22, 
and go = 0.6. Here again the agreement between theory and experiment 
is quite satisfactory except with regard to the nonlinear coefficient C3. For 
more quantitative comparisons the calculations need to incorporate rigid- 
rigid boundaries and the experimental geometry needs to be changed to ac- 
comodate many more periods, thereby better approximating the large-system 
limit. From the experimental values one obtains I + C1C3 > (the Newell cri- 
terion) indicating that the uniform traveling wave pattern is Benjamin- Feir 
stable. The stability of the traveling wave is bounded by the Eckhaus stable 
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wavenumber band 



1 + ClC 3 

Outside this stability boundary, the pattern is unstable to long wavelength 
perturbations and will undergo a subcritical bifurcation to a new pattern 
with a wavenumber within the stability boundary Given the experimentally 
obtained values for c\ and C3, we find that e# = 3£q(/c — k c ) 2 = 3em- In 
Fig. f|, we plot this Eckhaus boundary with finite size corrections [15]. This 



boundary agrees well with the experimental evidence that a transition took 
place around e = 0.0283 when gm = 0.0107. Since the transition was from 
m = 13tom = 12, both the amplitude and the precession frequency con- 
tained progressively more m — 12 component, thus becoming larger than 
if the state was purely mode 13, see Fig. |5|. We observed this variation 
during the transition. Since our cell is very uniform and the quasi one di- 
mensional wall convection structure is strictly periodic, there is virtually no 
long wavelength perturbation to force the instability to occur. This gives rise 
to long-lived "metastable" states (Zhong et at) when certain wall states have 
non-critical mode numbers that are outside the Eckhaus stability boundary. 



3.2 Bulk convection state 

As AT increased above the first bifurcation to the wall state, we observed a 
second transition in the Nusselt number measurements. Accompanied with 
this rise in heat transport, we also observed convective motion throughout 
the cell (the "bulk" state) and an increase in the thermal noise in the bottom 
plate which was supplied with a constant heat power. This indicates that 
the second transition was to time dependent convection, which was confirmed 
by analyzing time sequences of many shadowgraph images of the convection 
pattern. 

To characterize the bulk convection state, we measured its critical Rayleigh 
number R c b and critical wavenumber k c b for Q < 2145. The critical Rayleigh 
numbers, determined from heat transport measurements, are shown in Fig- 
ure [l]. Also plotted as a solid line is Chandrasekhar's linear stability analysis 
result for a laterally infinite system []I]] . The critical wavenumbers obtained 
from recorded shadowgraph images are shown in Figure || They also agree 
well with the results of Chandrasekhar, shown as the solid line. From each 
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Figure 6: The critical bulk wavenumber k c ^ vs f2 for the bulk state obtained 
from shadowgraph images. The solid line is from Chandrasekhar]!]]. 

shadowgraph image, we extracted their central portions, applied 2-D Fourier 
transformation with band pass frequency filtering and summed up the total 
power. Just above the onset where the shadowgraph image intensity variation 
is almost linear in the temperature gradient, the total power in the spectrum 
is J2i I A| 2 if the convection structures are composed of rolls of Aiexpiiki ■ x). 
In Figure [7] (a), (b) and (c), we show the background-divided and digitally- 
enhanced images of the convection patterns with Q = 2145 for values of R 
just below, nearly at and slightly above R c b, which was determined from the 
change in slope of the Nusselt number. The e values in the figure are de- 
fined relative to R c b. To more easily compare the bulk amplitude obtained 
from spectral analysis and heat transport data, we define an effective Nusselt 
number as the ratio of heat conducted by both convection states (bulk and 
side-wall) divided by the heat conducted by the side-wall state alone. Figure 
(d) shows this effective Nusselt number and J2i \Ai\ 2 versus e for Q = 2145. 
The linear relation between the total squared amplitude and e confirms that 
the bifurcation is forward. We conclude that the bulk convection state is the 
one given by Chandrasekhar 's linear stability analysis]!]]; because there is no 
appreciable amplitude of the side-wall state in the central portion the bulk 
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Figure 7: Shadowgraph images of convection pattern (a) below (R = 4.87 x 
10 5 , e = -0.062), (b) nearly at (R = 5.11 timeslO 5 , e = -0.016), and (c) 
above (R = 5.38 x 10 5 , e = 0.037) the onset of bulk convection, and (d) the 
effective Nusselt number Q/Q (•), ie. the ratio of the total heat conducted 
in the convecting state divided by the heat conducted by the side-wall state 
alone, and the squared amplitude (q) vs e. All the data are for Q = 2145 
and the e values for the patterns are labeled in (d). 
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state can grow from small amplitude despite being a second transition. 

Figure [I] shows that the bifurcation lines for the side wall state and the 
bulk convection state cross each other when 60 < Q < 70 for T = 2.5. For Q 
less than the intercept value, bulk convection sets in first and wall convection 
begins later. The bulk and wall convections are spatially separated, inter- 
acting only weakly. The cross-over point for the bifurcation lines is therefore 
not a co-dimension two point. A further point about this regime at small 
Q is that the small aspect ratio of the container may become important. 
Qualitatively nothing changes but for the quantitative comparison with the- 
ory these finite size effects are probably significant. In particular, the bulk 
state is the initial bifurcation and its presence might strongly influence the 
side-wall bifurcation which now happens at higher R; the high value, relative 
to theory, of the side-wall onset for fl = 16.93 suggests such an effect. The 
onset of the bulk state is also affected by the small aspect ratio when Q is 
small as evidenced by the small upward shift of R^ consistent with finite-size 
effects in non-rotating convection [|19|, |20| . Both of these shifts are visible in 
Fig □ 

There are several additional interesting aspects of the results obtained 
from shadowgraph images. First, there is a range of R around R c b where 
the convection pattern in the middle region consists of concentric rings. As 
is known|18[, concentric ring patterns in Rayleigh-Benard convection can 
be induced by thermal forcing from the side wall owing to the mismatch 
of thermal conductivities between the side wall and working fluid. In our 
case, the effective side wall for the bulk convection is the wall state, which is 
convecting rather vigorously at the second bifurcation point. For Q = 1090 
and near the onset of bulk convection (R/R cs = 2.246), Nu = 1.54. The 
active convection area is at most 50% of the cell, implying an effective thermal 
conductivity of the wall state at least twice that of the conducting water in 
the middle. The concentric rings eventually break up when e is increased 
and give rise to cellular structures. 

The total squared amplitude (J2i shown in Figure [7| (d), has a linear 

dependence on R. However, the extrapolation of this linearity to the base 
line yields a Rayleigh number slightly less than R c b, typically on the order of 
0.2%. Judging from the rounded corner of the Nu versus e curve near e = €2, 
this shift is caused by the appearance of the concentric ring pattern. 

A time sequence of the images taken at different points in parameter space 
reveals the source of the time dependence in the bulk convection state. For 
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as high as 2145, we observed a Kiippers-Lortz roll-orientation switching 
and/or turbine like rotation of convection cells. The dynamics stemming from 
the Kiippers-Lortz instability are aperiodic, which gives rise to the increased 
thermal noise in the bottom plate. This aspect of the result will be presented 
in greater detail elsewhere f2~T| . 



4 Conclusion 

In this paper, we presented experimental results for the initial bifurcations in 
rotating Rayleigh-Benard convection in small aspect ratio cylindrical cells. 
We confirmed a forward bifurcation from conduction to convection near rigid 
walls. Our experimental data agree beautifully with the theoretical results 
for the rotation dependence of the critical traveling wave frequency and mode 
number. We confirmed the prediction that the critical frequency goes linearly 
to zero as Q approaches zero. We also conjecture from the experiment that 
for insulating walls, the critical mode number will reach an asymptotic finite 
value as Q increases. We have shown for the second bifurcation that the 
convection states are very similar to the ones that are linearly unstable in 
the infinite system. Because of their spatial separation, the wall and bulk 
convection states interact only weakly in a narrow annular region between the 
two states. As aspect ratio increases, the separation becomes more distinct 
and the wall state becomes relatively more confined to the side wall. Since 
the wall state is quasi one dimensional, its contribution to the dynamics 
and transport properties of rotating thermal convection should diminish in 
the limit of large aspect ratio. The presence of such "surface" modes will, 
however, subtly affect aspects of the bulk nonlinear dynamics, as is evident 
by its forcing of the concentric ring pattern near the onset of bulk convection. 
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